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A Bayesian method application to the deconvolution ol EXAFS spectra is considered. It is shown 
that for purposes of EXAFS spectroscopy, from the infinitely large number of Bayesian solutions it 
is possible to determine an optimal range of solutions, any one from which is appropriate. Since 
this removes the requirement for the uniqueness of solution, it becomes possible to exclude the 
instrumental broadening and the lifetime broadening from EXAFS spectra. In addition, we propose 
several approaches to the determination of optimal Bayesian regularization parameter. The Bayesian 
deconvolution is compared with the deconvolution which uses the Fourier transform and optimal 
Wiener filtering. It is shown that XPS spectra could be in principle used for extraction of a one- 
electron absorptance. The amplitude correction factors obtained after deconvolution are considered 
and discussed. 
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I. INTRODUCTION 

The chief task of the extended x-ray absorption fine- 
structure (EXAFS) spectroscopy, determination of inter- 
atomic distances, rms fluctuations in bond lengths etc., 
is solved mainly by means of the fitting of parameterized 
theoretical curves to experimental ones. However, there 
exist obstacles for such a direct comparison: theory lim- 
itations and systematic errors. Among latter are various 
broadening effects. Fist of all, (i) this is the broadening 
arising from the finite energy selectivity of monochroma- 
tor and the finite angular size of the x-ray beam, (ii) 
The absorption even of strictly monochromatic x-ray ir- 
radiation by the electrons of a deep atomic level gives 
rise to photoelectrons with the finite energy dispersion 
due to the finite natural width of this level and the fi- 
nite lifetime of the core-hole, (iii) For x-ray energies far 
above the absorption edge the process of photoelectron 
creation (the outgoing from an absorbing atom) and the 
process of its propagation occur for essentially different 
time intervals. In other words, just created, the photo- 
electron 'does not know' where and how it will decay. 
Therefore the photoionization from the chosen atomic 
level and excitation of the remaining system can be con- 
sidered as independent processes, and hence the total ab- 
sorption cross-section, as a probability density of two in- 
dependent random processes, is given by the convolution 
of a one-electron cross-section and excitation spectrum 
W(AE). The latter is the probability density of the en- 
ergy AE capture at the electron-hole pair creation and 
is the quantity measured in x-ray photocmission spec- 
troscopy (XPS). For light elements there are examples of 
such deep and lengthy enough XPS spectra (see figure [l], 
taken from Q|). 

For all three cases the measured absorption coefficient 
jj, m is given by the convolution: 
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FIG. 1. From §: XPS spectra of MgO and MgF 2 in the 
vicinity of the Mg Is peak. Two secondary structures, due to 
plasmon losses and to double-electrton KL-zo, excitations, are 
detected. Zero of energy is placed at Mg K level (~ 1300 eV). 



where the broadening profile W(AE) and the meaning 
of the function fi depend on the considered problem. 
These can be, correspondingly: x-ray spectral density 
after the monochromator and the cross-section of ide- 
ally monochromatic irradiation; the Lorentzian function 
and the cross-section with a stationary initial level (of 
zero width); the excitation spectrum and a one-electron 
cross-section. 

It is the common practice in modern EXAFS spec- 
troscopy to account for the broadening processes (i)-(iii) 
at the stage of theoretical calculations by introducing 
into the one-electron scattering potential the imaginary 
correlation part which represents the average interaction 
between photoelectron and hole and their own polariza- 
tion cloud; in doing so, the choice of the correlation part 
is dictated by empiric considerations and can be different 
for different systems. 
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Another approach to the account for the broadening 
processes is to solve the integral equation (0) for the un- 
known fi. To find some solution of this equation is quite 
not difficult to do, the simplest way is to use the theorem 
about the Fourier transform (FT) of convolution. How- 
ever, it is known that the problem of deconvolution is 
an ill-posed one: it has an unstable solution or, in other 
words, the infinitely large number of solutions specified 
by different realizations of the noise. Thus, there is evi- 
dent necessity for determination of an optimal, in some 
sense, solution. Yet a less evident approach exists: to find 
an appropriate functional of the solution which itself be 
stable. 

A number of works have been addressed the problem 
of deconvolution, among them those concerning the x- 
ray absorption spectra. Loeffen et al. applied decon- 
volution with the Lorentzian function partly eliminating 
the core-hole life-time broadening. They used fast FT 
and the Wiener filter which is determined from the noise 
level which, in turn, is specified by the choice of the lim- 
iting FT frequency above which the signal is supposed 
to be less than noise. The arbitrariness of such a choice 
gives rise to rather different deconvolved spectra, which 
although remained obscured in [^j. 

Recently, for the deconvolution problem with a 
Lorentzian function Filipponi jj| used the FT and pro- 
posed the idea of the decomposition of an experimental 
spectrum into the sum of linear contribution, a special 
analytic function representing the edge and oscillating 
part. For the Lorentzian response, the deconvolution for 
the first two contribution is found analytically, for the 
latter one, numerically. The advance of such a decom- 
position is in that fact that now the FT of the oscillat- 
ing part is not dominated by the very strong signal of 
low frequency, therefore the combination of forward and 
backward FT gives less numerical errors. Notice that 
this method is solely suitable for the analytically given 
response. In addition, in || the choice of the filter func- 
tion (Gaussian curve) and its parameterization remained 
vague. Therefore the issue on the uniqueness or optimal- 
ity of the found solution was left open. 

In early work Q , for the solution of ill-posed problems 
the statistical approach was proposed. Following that 
work, in the present paper we shall consider the deconvo- 
lution problem in the framework of Bayesian method, de- 
tailed formalism of which was described in j^]. Since the 
parameterization is naturally involved in the Bayesian 
method, there exists a principle possibility to choose an 
optimal, in some sense, parameter. Here we shall scru- 
tinize the problem of such a choice which is relevant to 
any spectroscopy. We shall show that this problem is ab- 
sent in EXAFS spectroscopy because though the EXAFS 
spectrum itself does depend on the regularization param- 
eter, its FT does not in the range of real space used for the 
analysis. In Sec. || we discuss the choice of the optimal 
deconvolution for a model Gaussian response and com- 
pare the results of Bayesian approach with the results of 
FT combined with Wiener filtering. In Sec. Ill we utilize 



the deconvolution to an experimental spectrum in order 
to eliminate the aforementioned broadening processes. 



II. THE CHOICE OF OPTIMAL 
DECONVOLUTION 

First of all, we shall show that the deconvolution prob- 
lem really has the infinitely large number of solutions. 
In the present paper we use for examples the absorp- 
tion spectrum of Ndi.ssCeo.isCuO^a above Cu K edge 
collected at 8 K in transmission mode at LURE (beam- 
line D-21) using Si(lll) monochromator and harmonics 
rejecting mirror; energy step ~ 2 eV, total amount of 
points 826 (from 8850 eV to 10500 eV), each one recorded 
with integration time of 10 s. Let us take for a while for 
the response function a simple model form: W(E) = 
Cexp(— E 2 /2T 2 ), where C normalizes W to unity, T is 
chosen to be equal to 4 eV. 

In [El we showed how to construct a regularized solu- 
tion of the convolution equation in the framework of the 
Bayesian approach. For that, one needs to find eigenval- 
ues and eigenvectors of a special symmetric N x N matrix 
determined by the experimental spectrum, N is the num- 
ber of experimental points. Using that approach, find a 
solution fJ,(E) for an arbitrary regularization parameter 
a and perform its convolution with W . The /t m obtained, 
ideally, must coincide with /x m . Introduce the character- 
istics of the solution quality, the normalized difference of 
these curves: 
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regularizer a 

FIG. 2. The quality of the deconvolution R 2 vs. regulariza- 
tion parameter a (solid) . Dashed lines — signal-to-noise ratio 
before and after deconvolution. Two peaks — the posterior 
density functions p(ajd) (left) and p(a|d, a) (right) for the 
regularization parameter a. Straight line — the determinant 
of Bayesian matrix as a function of a. 



2 



(2) 



16 20 



where the summation is done over all experimental 
points. Figure || shows the dependence R 2 on a. That 
fact that the quality of the found solutions is practically 
the same for all a < 1 is a clear manifestation of ill— 
posedness of the problem: there is no a unique solution. 
How to chose an optimal one? It turns out, that for pur- 
poses of EXAFS spectroscopy there is no need of that and 
any solution from the optimal range (here, a < 1) is suit- 
able. At arbitrary a from the optimal range found the de- 
convolution n(E), extract the EXAFS function x(fc) • k w 
in a conventional way, where k is the photoelectron wave 
number, and find its FT. In figure ||we show the EXAFS- 
functions obtained after the Bayesian deconvolution with 
a = 1 and a = 0.01, and their FT's. As seen, although 
the EXAFS-function itself does depend on a, its FT prac- 
tically does not. Thus, if one uses for fitting a range of 
r-space (in our example, up to r max = 8 A) or filtered k- 
space, the problem of search for the optimal a is not rele- 
vant. Nevertheless, below we propose several approaches 
to the solution of this problem, for instance for XANES 
spectroscopy needs. 

(1) For the regularization parameter a itself one can 
introduce the posterior probability density function [|]]^] 
and choose a with a maximum probability density. It can 
be done either using the most probable value of noise 
or assuming the standard deviation a of the noise of 
the absorption coefficient to be known (for our spectrum 
a = 9 • 10~ 4 , as determined from the FT following ||). 
In figure these probability densities are drawn as, cor- 
respondingly, p(a|d) and p(a|d, a), and their most prob- 
able values are a\ = 0.021 and a\ a = 0.044. 

(2) The optimal regularization can be determined 
from the consideration of the signal-to- noise ratio S/N . 
The Shannon-Hartley theorem states that / max = 
Bin (1 + S/N), where / max is the maximum informa- 
tion rate, B is the bandwidth. The authors of Q are of 
opinion that deconvolution is a mathematical operation 
that conserves information. Therefore from the theorem 
follows that one pays for an increase in bandwidth, re- 
sulted from deconvolution, via a reduction in S/N ratio. 
The thesis on 7 max conservation is quite questionable, 
since for deconvolution one should introduce additional 
independent information about the profile of broaden- 
ing. What quantity is conserved in deconvolution is hard 
to tell. Here to the contrary, for the optimal a we de- 
mand to conserve S/N . Define S/N as the ratio of mean 
values of the EXAFS power spectrum over two regions, 
r < 15 A and 15 < r < 25 A. The regularization param- 
eter at which the S/N is conserved is denoted in figure 
U as a?2 = 5.54. The signal-to-noise ratio can be defined 
in a different way. Since the Bayesian methods work 
in terms of posterior density functions, for each experi- 
mental point one can find not only the mean deconvolved 
value but also the standard deviation 5/idcconv from which 
one finds 8\ = ^deconv/Vcij where no is the atomic-like 




FIG. 3. % ' ^ 2 obtained without deconvolution and after 
that with a = 1 and a — 0.01. In the middle — the envelope 
of the initial x an d rms deviations of the deconvolved values 
(dots). Below — the absolute values of the Fourier transform 
(the dashed and the dotted curves practically merge). 



absorption coefficient constructed at the stage of EXAFS 
function extraction. It is reasonable to compare 8\ values 
with the envelope of EXAFS spectrum (figure ^ middle) . 
As seen, at small a the noise dominates over the signal 
in the extended part of the spectrum. The regularization 
parameter at which they match is the optimal one, a.2- 

(3) For the Bayesian deconvolution it is necessary to 
find eigenvalues and eigenvectors of a special symmetric 
matrix g. It turns out that the determinant of this matrix 
varies with a over hundreds orders of magnitude. At 
small a's the matrix is poorly defined, large a's yield very 
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FIG. 4. Module of the FT of initial /x m and deconvolved fi 
at different a. Below — filters transforming /im(i~) to /Lt(r) 
obtained after Bayesian deconvolution (thin solid lines), after 
deconvolution based on the FT (dashed line), and after de- 
convolution based on the combination of the FT and Wiener 
filtering (thick solid line). 



large det(g) (figure ||). Both cases give large numerical 
errors because of ratios of very small or very large values 
in calculations. As an optimal parameter we choose = 
1.41 at which det(g) ~ 1. 

The cases (1) and (2) require to determine the noise 
level, which demands additional variables (for instance, 
the limiting frequencies of FT). The case (3) does not 
explicitly concern the noise. Due to the dependence of 
lg[det((7)] on a appears to be linear, which readily allows 
one to find the optimal parameter, the case (3) is more 
preferable from the practical point of view. Below, for 
deconvolution of the real broadening processes we use 
the optimal parameter a^. 

It is of certain interest to compare the Bayesian method 
of deconvolution with the widely known method combin- 
ing optimal Wiener filtering and the convolution theo- 
rem, where the conjugate variables of the FT are not k 
and 2r adopted in EXAFS but E and r. According to 
the theorem, yU m (r) = W(t) ■ /i(r), where for our model 
Gaussian response W(t) = cxp(r 2 r 2 /2). A simple back 
FT of the ratio fi m (T)/W(r) will give the thought solu- 
tion n(E) but very noisy. Therefore fJ-mir) at large t has 
to be smoothed. Figure f| shows module of the FT of 
the measured spectrum and of the Bayesian deconvolved 
spectra. The latter are merged at r < 0.25 eV -1 . In the 
bottom part of the figure the ratios |FT/Udeconv|/|FT^ m | 



are shown for different a. As seen, the Bayesian decon- 
volution performs the effective filtration of spectra with 
the limiting frequency r max depending on a. 

The optimal, in the least-square sense, Wiener filter 
is expressed as §: $(t) = (1 + \n(r)\ 2 /IFT^] 2 )- 1 , 
where |n(r)| 2 is the power spectrum of the noise re- 
placed here by 0.01, the mean value of |m(t)| 2 over the 
range r > 0.4 eV -1 . As seen in figure [|, the effective 
Wiener filter W(t)& transforming ^m(r) to ^(r) is close 
to the effective filter of the Bayesian deconvolution with 
a = 1 ps a^. Notice, however, that the limiting fre- 
quency for the estimation of noise power spectrum was 
chosen rather arbitrarily. 

In closing this section, it should be noticed that apart 
from the possibility of determination of the deconvolu- 
tion errors and the possibility of the optimal regulariza- 
tion parameter choice, the Bayesian deconvolution has 
the advantage of the capability to take into account a pri- 
ori information about the smoothness and shape of the 
solution (see details in ||). In addition, in the Bayesian 
method the response function W(E — E') could be of 
more general form W(E — E',E), which will be useful 
for deconvolution of the instrumental broadening because 
monochromator energy resolution depends noticeably on 
the angular position and, hence, on the energy of the 
output x-ray beam. 



III. APPLICATIONS OF DECONVOLUTION 

We have seen that the Bayesian method proves to be 
effective for deconvolution of EXAFS spectra, and the 
choice of the regularization parameter appears to be ir- 
relevant. Now we perform the deconvolution of various 
types of broadening, for which purpose specify the corre- 
sponding response functions W(E — E' ,E). 



A. Instrumental broadening 

The monochromator resolution is determined by the 
rocking curve width 59b and by the vertical beam diver- 
gence 59 j_. For the monochromator Si(lll) at E = 9keV 
the rocking curve width is 89 B = 32.4 ^irad (FWHM) 
§, the beam divergence (LURE, D-21) <56»_l = 150 ^rad. 
Strictly speaking, the resulting spectral distribution is 
given by the convolution of rocking curve and the an- 
gular beam profile. But since 89 b -C 56±, the energy 
selectivity is determined by 89 ±, namely: 8E±/E = 
59 x cot 9 B = 59 ± ^/(2Ed/ch) 2 - 1, where 9 B is Bragg 
angle, d is Bragg plane spacing. Modelling the spectral 
distribution by a Gaussian function, one obtains: 



W(E — E', E) cxexp 



(E - E' f 



2*1(25) 



,<r±(E) = 



5Ej_(E) 
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where the normalization constant must be calculated at 
each E value. For our sample spectrum, uj_(8850eV) = 
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2.46 eV and a± (10500 eV) = 3.49 eV. 



B. Lifetime broadening 

For deconvolution of the lifetime broadening described 
by a Lorentzian function W(AE) oc [(AE/T K ) 2 + l]" 1 , 
we take as the initial spectrum /x m the spectrum /ii ns t r ob- 
tained after the instrumental deconvolution. According 
to f§, the width of Cu K level (FWHM) equals 1.55 eV, 
from where T K = 0.775 eV. 



C. Multielectron broadening 

There are certain difficulties in measuring XPS spec- 
tra near (and deeper) the deepest atomic levels: the 
monochromatic x-ray sources of high energy are required; 
for long enough spectra (~ 100 eV) a photoelectron an- 
alyzer with a broad energy window and long integration 
time are necessary. Unfortunately, for lack of experi- 
mental XPS spectra in cuprates near Cu Is level, we 
can use a model representation of the response W(AE). 
For that we take the estimations of position, inten sity 
and width of the secondary KM23 excitation from p0[ : 



Ekm 23 — Ek = 85 eV; Ikm 23 /Ik = 0.03; Tkm 2 
finally, for the broadening function we have: 
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Again, as the initial spectrum we take the spectrum 
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FIG. 5. Module of the FT of various EXAFS functions: 
initial; obtained from the instrumentally deconvolved absorp- 
tance; the latter was deconvolved with the Lorentzian re- 
sponse and with the total multielectron response, these two 
give the EXAFS FT's practically merged. Bottom: the ratio 
of amplitudes of initial x an d deconvolved x- The Sq value 
was calculated from atomic overlap integrals. 



IV. DISCUSSION 

With the specified response functions W, perform the 
Bayesian deconvolution of the absorption coefficient at 
the optimal regularization parameter, then construct the 
EXAFS function for which calculate FT and the ampli- 
tudes and phases (see figure ||). Just as for the model re- 
sponse in Sec. [il], the deconvolution leads to the increase 
of EXAFS oscillations. As appeared, the deconvolution 
has practically no influence on the first FT peak origi- 
nating from the shortest scattering path. It is clear why: 
the oscillations corresponding to this peak are essentially 
wider then the response W (for these, W is almost a 8- 
function), and this is more true for the extended part of 
a spectrum, due to the period of the oscillations there 
is even longer (in _E-space). Thus, it is in the extended 
part where (M, the solution of equation ([[]), less differs 
from fi m . 

In modern EXAFS spectroscopy the difference between 
amplitudes of experimental and calculated spectra are 
taken into account by the reduction factor Sq which is 
either treated as a fitting parameter or estimated from 



the relaxation of the core-hole as the many-electron over- 
lap integral. In many works this factor is considered to 
be independent from energy, however, as noted in review 
by Rehr and Albers |11|, it must be path-dependent and 
energy-dependent. At the bottom of figure ^| we draw 
the ratios A(k) of the amplitude of initial EXAFS spec- 
trum to that of the deconvolved one. Here, they were 
calculated relatively both Xm and Xi ns tr ■ For comparison 
we show the factor Sq as computed by feff code [ ^2[ . 
At large fc, where the noise become comparable with the 
EXAFS signal, the ratios A(k) have significant errors. 
However, the general trend of the curves corresponds to 
the expected one [ p3[ : A{k) is minimal at intermediate 
EXAFS energies, while at both low and high energy A(k) 
reduces to unity. In addition, there are some phase shifts 
between the initial spectrum and the deconvolved ones. 
But these shifts are found to be quite small: less than 
0.2 rad at k < 4 A and less than 0.1 rad at k > 4 A. 

The Lorentzian broadening of the EXAFS spectrum 
with a half-width T is similar to the effect of the imagi- 
nary part of the self-energy with Im S = T. The resulting 
reduction factors, i.e. the ratios the of amplitudes calcu- 
lated with and without the imaginary part, are analogous 
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to the reduction factors obtained by us relatively Xinstiv 
A(k) = am(xinstr)/am(xdcconv)- However, up to now the 
reduction factors relatively measured EXAFS spectrum 
were considered, A(k) = am(x m )/arn(xdeconv)- As seen 
(figure ||), these are the noticeably different quantities. 
That is why for the correct analysis and comparison of 
spectra taken at different experimental conditions, the 
instrumental deconvolution must be the first step. 

For our example spectrum and the chosen response 
functions, deconvolution of the lifetime broadening and 
deconvolution of the multielectron broadening are practi- 
cally undistinguishable (figure ||, top), i.e. the secondary 
weak peak in the excitation spectrum W(AE) has very 
little effect. The main effect of using real excitation 
spectra is expected from the presence and taking into 
account the plasmon losses which have a considerable 
integral weig ht (figure 0). Because of their very broad 
spectral distribution, their effect consists in the change 
of the EXAFS spectrum as a whole. In the present paper 
this contribution was not taken into account for lack of 
appropriate experimental information. 

Near the absorption edge, where the photoelectron ki- 
netic energy is low, the core-hole relaxation processes are 
of certain importance for the photoelectron propagation. 
Here we do not consider the validity of the neglect of this 
effect, but refer to the review |pT[. 



V. CONCLUSION 

To take into account the many-electron effects, there 
exist, in principle, two approaches: (a) to include into 
a one-electron theory relevant amendments or (b) to ex- 
tract a one-electron absorptance from the total one and 
to use then a pure one-electron theory. The first, tra- 
ditional, approach invokes semi-empirical rules, but not 
ab initio calculations, to construct the exchange corre- 
lation part of the scattering potential, with the empiri- 
cism being based on the comparison with experimental 
spectra already broadened. In the present paper we have 
shown the principle way for the second approach, using 
the solution of integral convolution equation, the kernel 
in which is the excitation spectrum measured in XPS 
spectroscopy. Notice, that owing to the specific way of 
the structural information extraction from the EXAFS 
spectra, in which an isolated signal in r-space or a fil- 
tered signal in &;-space is used, we have not committed 
a sin against the fact that the integral convolution equa- 
tion is an ill-posed problem, because from the infinitely 
large number of solutions it is possible to determine an 
optimal range, also infinitely large, of solutions any one 
from which is appropriate. 

Because of some technical difficulties, it is impossible 
so far to measure XPS spectra near deep core levels. 
Therefore, the desirable pure one-electron absorptance 
is unavailable. Nevertheless, as we have shown, it is pos- 
sible to perform an accurate instrumental deconvolution 



and deconvolution of the lifetime broadening. These pro- 
cedures make the comparison between calculated and ex- 
perimental spectra more immediate and the final results 
of EXAFS spectroscopy more reliable. 

All the stages of EXAFS spectra processing including 
those described here are realized in the freeware program 
viper ITl. 



ACKNOWLEDGMENTS 

The example spectrum was measured by Prof. A. P. 
Mcnushcnkov. The author wishes to thank Dr. A. V. 
Kuznetsov for many valuable comments and advices. The 
work was supported in part by RFBR grant No. 99-02- 
17343. 



[1] A. Di Cicco, M. De Crescenzi, R. Bernardini, and G. 

Mancini, Phys. Rev. B 49, 2226 (1994). 
[2] P. W. Loeffen, R. F. Pettifer, S. Miillender, M. A. van 

Veenendaal, J. Roler, and D. S. Sivia, Phys. Rev. B 54, 

14877 (1996). 

[3] A. Filipponi, J. Phys. B: At. Mol. Opt. Phys. 33, 2835 

(2000) . 

[4] V. F. Turchin, V. P. Kozlov, and M. S. Malkevich, Sov. 

Phys. Usp. 13, 681 (1971). 
[5] K. V. Klementev, J. Phys. D: Appl. Phys. 34, accepted 

(2001) ; e-arXiv:physics/0003086. 

[6] M. Newville, B. I. Boyanov, and D. E. Sayers, J. Syn- 
chrotron Rad. 6, 264 (1999), (Proc. of Int. Conf. XAFS 
X). 

[7] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. 

Flannery, Numerical Recipes in Fortran 77: The Art of 

Scientific Computing, Second Edition (Cambridge Univ. 

Press, Cambridge, 1992), chap. 13.3. 
[8] M. Sanchez del Ri'o, C. Ferrero, and V. Mocella, SPIE 

proceedings 3151, 312 (1997). 
[9] M. O. Krause and J. H. Oliver, J. Phys. Chem. Ref. Data 

8, 329 (1979). 

[10] A. Di Cicco and F. Sperandini, Physica C 258, 349 
(1996). 

[11] J. J. Rehr and R. C. Albers, Rev. Mod. Phys. 72, 621 
(2000). 

[12] A. L. Ankudinov, B. Ravel, J. J. Rehr, and S. D. Con- 

radson, Phys. Rev. B 58, 7565 (1998). 
[13] J. J. Rehr, E. A Stern, R. L. Martin, and E. R. Davidson, 

Phys. Rev. B 17, 560 (1978). 
[14] K. V. Klementev, 2000, VIPER for Windows (Visual 

Processing in EXAFS Researches), 

freeware, www.crosswinds.net/~klmn/viper.html . 



G 



